Estimates of the uncertainties associated with models 
of the nucleon structure functions in the A production region 
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Theoretical studies of the inclusive electron-nucleus cross section at beam energies up to few GeV 
show that, while the region of the quasi-elastic peak is understood at quantitative level, the data 
in the A production region are sizably underestimated. We analize the uncertainty associated with 
the description of the nucleon structure functions W\ and Wi and its impact on the nuclear cross 
section. The results of our study suggest that the failure to reproduce the data is to be mostly 
ascribed to the poor knowledge of the neutron structure functions at low Q 2 . 
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In view of the rapid development of neutrino physics, 
leading to significant improvements in the experimental 
accuracy, the treatment of nuclear effects in data analysis 
is now regarded as one of the main sources of systematic 
uncertainty 0,0|- 

Much of the information needed to understand the 
nuclear response to neutrino interactions in the energy 
range E v — 0.5—3 GeV, relevant to many neutrino exper- 
iments, can be extracted from the results of experimental 
and theoretical studies of electron-nucleus scattering (for 
a recent review see, e.g., Ref. |3j). In this kinematical 
regime, in which single nucleon knock out is known to be 
the dominant reaction mechanism, both quasi-elastic and 
inelastic processes, leading to the appearance of hadrons 
other than protons and neutrons, must be taken into ac- 
count. 

Non relativistic nuclear many-body theory (NMBT) 
provides a consistent framework, suitable to describe 
electron-nucleus interactions in a variety of kinematical 
conditions, ranging from quasi-elastic to deep inelastic 
scattering (see, e.g., Ref. |i(). In the impulse approx- 
imation regime, in which scattering off a nuclear tar- 
get reduces to the incoherent sum of elementary pro- 
cesses involving individual nucleons, the basic elements 
of the NMBT approach are the nucleon spectral function 
P(p,E), yielding the energy and momentum distribu- 
tion of the knocked out nucleon, and the electron-nucleon 
cross section, that can be written in terms of the two 
structure functions W^ and W^ (N —p,n). 

A study of inclusive electron scattering off oxygen at 
beam energies between 0.7 and 1.2 GeV, carried out 
within NMBT J5J], has recently shown that, while the 
data in the region of the quasielastic peak are accounted 
for with an accuracy better than ~ 10 %, theory fails to 
explain the measured cross sections at larger electron en- 
ergy loss, where A production dominates. Based on the 
fact that the calculation of the cross section at the quasi- 
elastic and A production peak involves integrations of 
P(p, E) extending over regions of the (p, E) plane which 
almost exactly overlap one another, the authors of Ref. 



|5J argued that the disagreement between theory and data 
is unlikely to be imputable to deficiencies of the spectral 
function, and should rather be ascribed to the description 
of the elementary electron-nucleon cross section. 

This short note is aimed at providing a quantitative 
estimate of the uncertainty associated with the available 
models of the nucleon structure functions at low Q 2 . The 
impact on the nuclear cross section is also analyzed com- 
paring theoretical results to the data of Refs. @, 0- 

The calculations of Ref. [5j were carried out using 
the Hohler-Brash parameterization of the nucleon form 
factors [3, 13] , resulting from a fit which includes the re- 
cent Jefferson Lab data [l(| , and the Bodek and Ritchie 
(BR) parametrization of the proton and neutron inelastic 
structure functions pd| , covering both the resonance and 
deep inelastic region. 

The structure functions of Ref. [IlJ have been obtained 
from a global fit to the electron-proton and electron- 
deuteron cross sections measured at SLAC [12|, span- 
ning the kinematical domain 1 < Q 2 < 20 GeV 2 and 
0.1 ;$ x < .75, x being the Bjorkcn scaling variable. As a 
consequence, as pointed out in Ref. [5j , using the results 
of Bodek and Ritchie in the kinematical region covered by 
the data of Ref. Q, corresponding to Q 2 < 0.2 GeV 2 at 
the A production peak, involves an extrapolation whose 
validity needs to be carefully investigated. 

Additional uncertainty is associated with the neutron 
structure functions Wf and W% , that have been ex- 
tracted from deuteron data subtracting the proton contri- 
bution and unfolding nuclear effects. The authors of Ref. 
[ill ] followed the approach of Atwood and West [lj| , with 
the deuteron wave function obtained using the Hamada- 
Jonston (HJ) model of the nucleon-nucleon (NN) poten- 
tial Q. 

To gauge the systematic error involved in the determi- 
nation of Wi and W^ > we have compared the electron- 
deuteron cross sections computed using the HJ deuteron 
wave function to those obtained using the state-of-the- 
art NN potential of Ref. [15|. generally referred to as 
Argonne v±s (A18) potential. The HJ and A18 models 



significantly differ in the description of the short range 
repulsive component of the NN force. The HJ exhibits a 
hard core of radius r c ~ 0.5 fm, leading to a vanishing 
deuteron wave function at r < r c , while the A18 poten- 
tial smoothly reaches a large but finite positive value at 
r = 0. The differences in the short range behavior of 
the wave functions in coordinate space are reflected by 
sizable differences in the momentum distributions n(p) 
at |p| £ 1 GeV. However, as n(p) falls by a factor ~ 10 6 
in the range < |p| < 1 GeV, these differences have 
negligible effect on the calculated electron-deuteron cross 
section. For example, at beam energy E e = 4 GeV and 
electron scattering angle 8 e = 30° the results obtained 
using the HJ and A18 momentum distributions differ by 
at most 2 % over the energy loss range extending from 
pion production threshold, corresponding to v ~ 1.05 
GeV, to v ~ 2.5 GeV. 

We have also investigated the effect of the ambiguity 
implied by the relativistic normalization of the momen- 
tum distribution, which is obtained from the nonrela- 
tivistic wave function in momentum space. A number of 
theoretical calculations have been carried out using the 
normalization of Ref. E3 instead of the one suggested 
by Attwood and West, which amounts to replacing 



n(p) 



vV 



n(p) , 



(1) 



m being the nucleon mass. Our results show that, as the 
momentum distribution is sharply peaked at small |p|, 
the presence of the additional normalization factor does 
not produce any appreciable effects. 

Having established that the uncertainty arising from 
the treatment of nuclear effects is small, we can now ad- 
dress the question of whether the fit of Ref. E3 provides 
a reasonable description of the structure functions at low 
Q 2 . To clarify this issue we have compared the calculated 
electron-deuteron cross section, obtained using the A18 
momentum distribution and the BR structure functions, 
to the Jefferson Lab data of Refs. [l7lE^ at E„ = 2.445 
GeV and 6 e — 20°. This kinematical setup corresponds 
to the lowest Q 2 at the A production peak available in 
the Jefferson Lab dataset, namely Q 2 = 0.54 GeV 2 . 

The results, represented by the solid line in the lower 
panel of Fig. Q show that in the A region the measured 
cross section is underestimated by as much as ~ 25 %. 
On the other hand, the results displayed in the upper 
panel show that the BR model provides a much better 
account of the electron-proton cross section, thus imply- 
ing that the disagreement with deuteron data is to be 
mostly ascribed to the neutron contribution. 

For comparison, in the upper panel we also include 
the cross sections calculated using the structure functions 
resulting from a recent fit including the Jlab data (19, 20] 
and from the dynamical model of Refs. |2lLl22| . It clearly 
appears that the different models describe the data in 
the region of the A peak with comparable accuracy. The 
distinctive behavior exhibited by the predictions of the 
model of Refs. |2l], [23 , at both lower and higher energy 
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FIG. 1: (Color online) Electron-proton (upper panel) and 
electron-deuteron (lower panel) cross sections at beam en- 
ergy 2.445 GeV and scattering angle 20°, corresponding to 
Q 2 = 0.54 GeV 2 at the A production peak. The solid lines 
show the results of calculations carried out using the structure 
fumctions of Ref. [llj , while the dashed and dot-dash lines in 
the upper panel have been obtained u sing the proton struc- 
ture functions of Refs. |19ll20| and |2ll l22j . respectively. The 
dashed lines labelled n and p in the lower panel correspond to 
the neutron and proton contributions to the deuteron cross 
section, respectively. The experimental data are taken from 
Refs. El Eg. 



loss, is to be ascribed to the fact that it does not include 
the contributions of other resonances and non resonant 
pion production. 

In Refs. [llL Il2j the deuteron cross section is written 
in the form 



(Td = Vp + er„ 



(2) 



where the proton and neutron contributions, er p „ can 
in turn be expressed in terms of the proton and neu- 
tron structure functions, smeared by nuclear binding and 
Fermi motion. Using the Atwood-West formalism and 
the measured electron-proton cross sections, a p , one can 
then compute the ratio 



Om 



(3) 



whose value provides a measure of nuclear effects. Note 
that Eq. JSJ) is based on the premise that the ratios 
between the unsmeared and smeared structure functions, 
Wf 2 and Wf 2 , satisfy E2 
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The results of numerical calculations carried out in the 
A production region show that S p , while varying sharply 
as a function of the invariant mass of the hadronic final 
state, W, exhibits a rather weak dependence on Q 2 . 




FIG. 2: (Color online) The nucleon structure functions W 2 N 
(N = n,p) at E e = 2.445 GeV and 9 e = 20°, corresponding 
to Q 2 = 0.54 GeV 2 at the A production peak, plotted as a 
function of the invariant mass of the hadronic final state. The 
shaded area represents the W% resulting from our analysis, 
while the solid and dashed line correspond to W% and W% of 
Ref. Flfl . respectively. 

Assuming that the smearing ratios for the proton and 
neutron cross section, S p and S n , be close to one an- 
other, the unsmeared neutron cross section can finally be 
estimated from 



S p (Jd 



(5) 



The neutron structure function obtained applying the 
above procedure at beam energy E e = 2.445 GeV and 
electron scattering angle 9 e = 20°, plotted as a function 
of IT, is shown by the shaded area in Fig. |2 The result- 
ing W£ ', which accounts for the measured deuteron cross 
sections by construction, turns out to be significantly 
larger than the one obtained from the fit of Ref. [llj | . 
the difference at the A production peak being ~ 60 %. 
For comparison, the proton structure function Wf is also 
shown. 

To gauge the impact of the uncertainty in the neutron 
structure functions on the nuclear cross section, we have 
analyzed electron scattering off carbon at E e — 1.3 GeV 
and 6 e = 37.5° and compared our results to the SLAC 
data of Ref. @|. 

The shaded area in the upper panel of Fig. [3] shows 
the results of our calculations, carried out within the 
Plane Wave Impulse Approximation (PWIA) formalism 
described in Ref. [5j with the spectral function obtained 
from the Local Density Approximation (LDA) [23. We 
have used the proton structure functions of Ref. [llj . 
while for the neutron the BR fit has been modified ac- 
cording to 



W%{Q\W) - S(W)W?(Q 2 ,W) , 



(6) 
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S(W) being the ratio between the results of the analysis 
described in the previous section (shaded area of Fig. 
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FIG. 3: (Color online) Upper panel: inclusive electron 
scattering cross section off carbon at E e = 1.3 GeV and 
6 e — 37.5°, corresponding to Q 2 = 0.4 GeV 2 at the A pro- 
duction peak, as a function of the electron energy loss v. The 
solid line corresponds to theoretical calculations carried out 
using the proton and neutron structure functions of Ref. [Till , 
while the shaded region has been obtained replacing the BR 
W\ and W2 with those resulting from our analysis. The data 
are taken from Ref. [g. Lower panel: same as in the upper 
panel, but for oxygen target, E e = 1.2 GeV and 9 e = 32°, 
corresponding to Q 2 — 0.26 GeV 2 at the A production peak. 
The data are taken from Ref. 01 ■ 



|2J) and the corresponding BR fit (solid line of Fig. |2J). 
In view of the fact that in the kinematical setup of the 
carbon data the value of Q 2 at the A production peak 
is - 0.4 GeV 2 , using S(W) extracted at Q 2 - 0.5 GeV 2 
appears to be reasonable. 

Comparison with the solid line, obtained using the BR 
fit for both the proton and the neutron, shows that the 
neutron structure functions that reproduce the deuteron 
data at Q 2 = 0.54 GeV 2 also provide a much better de- 
scription of carbon data at electron energy loss v > 0.6 
GeV. Discrepancies between theory and experiment still 
persist at lower v, in the region of the dip between the 
quasi-elastic and the A production peak. 

The inclusive cross section in the dip region is known 
to be affected by processes in which the electron couples 
to the two-body component of the nuclear electromag- 
netic current |24| . arising from 7r and p meson exchange 
between nucleons. However, the amount of strength asso- 
ciated with meson exchange currents (MEC) is still some- 
what controversial, and the possibility that the cross sec- 
tion in the dip region might be ascribed to an asymme- 
try in the shape of the quasi-elastic peak has also been 
suggested (see, e.g., J25J and references therein). The 
contribution of MEC is not taken into account in the 



PWIA approach of Ref. J5j, and its significance in the 
kinematical regime under discussion should be carefully 
investigated. 

The lower panel of Fig. [SJshows a comparison between 
our theoretical results and the oxygen data discussed in 
Ref. ya|- It clearly appears that at the low Q 2 of the 
data of Ref. Q the enhancement resulting from using the 
neutron structure functions obtained from our analysis is 
not sufficient to explain measured cross section. 

The results discussed in this paper show that the 
structure functions obtained from the global fit of Refs. 
|lltll2j . while providing a quantitative account of proton 
data, fail to explain the measured deuteron cross sections 
at low Q 2 (< 1 GeV 2 ). 

Employing a simple and somewhat crude procedure, 
similar to the one followed by Bodek and Ritchie, we 
have extracted the neutron structure functions at Q 2 ~ 
0.5 GeV 2 from the Jefferson Lab electron-proton and 
electron-deuteron data. The resulting W™ and W% turn 
out to be much larger than those of Ref. [ll| in the region 
of the A production peak. 

When used to calculate the nuclear cross section in 
the impulse approximation scheme, the neutron struc- 
ture functions obtained from our analysis provide a sat- 
isfactory description of carbon data at Q 2 = 0.4 GeV 2 . 
On the other hand, oxygen data at lower Q 2 are still 
sizably underestimated. This feature is likely to reflect 
an appreciable Q 2 -dependence of the function S(W), de- 
scribing the difference between our results and the BR 
fit. 

The results of Figs^ an d El have been obtained ne- 



glecting final state interactions (FSI), which are known 
to lead to a redistribution of the inclusive strength [j]. 
Within the approach of Ref. [n| , FSI effects in the quasi- 
elastic channel are included through a folding procedure 
that accounts for nucleon-nucleon rescattering processes. 
However, the extension to inelastic channels involves ad- 
ditional problems and has not been carried out yet. To 
estimate the relevance of FSI in the A production region 
we have folded the PWIA inelastic cross sections of Figs. 
UandElusing the folding functions obtained from the ap- 
proach of Ref. 16] . The results, showing that the main 
effect is a quenching of the peak of less than ~ 2% and 
4% in deuteron and carbon, respectively, suggest that 
FSI do not significantly affect our analysis. 

The ability of dynamical models [U IH |H |27| to 

explain the nuclear cross section should also be investi- 
gated. The results of a recent calculation [23 based on 
the model of Ref. Eg, are in fairly good agreement with 
the data of Ref. Q in the A production region at Q 2 ~ 
0.2 GeV 2 . However, the fact that measured quasi-elastic 
cross section is sizably overestimated seems to point to 
deficiencies in the treatment of nuclear effects. 

In conclusion, our study suggests that the neutron 
structure functions in the A production region at low 
Q 2 , which is still poorly known, may be extracted from 
a fit to the upcoming deuteron data from Jefferson Lab 
|28j . Valuable complementary information will also come 
from the direct measurement of the free neutron struc- 
ture functions at 1 < Q 2 < 5, presently being carried out 
by the Jefferson Lab E03-12 (BoNuS) collaboration J2<|. 



[1] Proceedings of NuIntOl, Eds. J.G. Morfin, M. Sakuda [17] 

and Y. Suzuki. Nucl. Phys. B (Proc. Suppl.) 112 (2002). 
[2] Proceedings of Nulnt04, Eds. F. Cavanna, P. Lipari, C. [18] 

Keppel and M. Sakuda. Nucl. Phys. B (Proc. Suppl.) 139 [19] 

(2005). 

[3] O. Benhar, D. Day and I. Sick, nucl-ex/0603029 Submit- [20] 
ted for publication in Reviews of Modern Physics. [21] 

[4] O. Benhar, V.R. Pandharipande and S.C. Pieper. Rev. 

Mod. Phys. 65 (1993) 817. [22] 

[5] O. Benhar, N. Farina, H. Nakamura, M. Sakuda and R. 

Seki, Phys. Rev. D 72 (2005) 053005. 
[6] R.M. Sealock et al, Phys. Rev. Lett. 62 (1989) 1350. [23] 

[7] M. Anghinolfi et al, Nucl. Phys. A 602 (1996) 405. 
[8] G. Hohler et al, Nucl. Phys. B114 (1976) 505. [24] 

[9] E.J. Brash, A. Kozlov, Sh. Li, and G.M. Huber, Phys. 

Rev. C 65, 051001(R) (2002) 051001 (R). [25] 

[10] M.K. Jones et al, Phys. Rev. Lett. 84 (2000) 1398. 
[11] A. Bodek and J.L. Ritchie, Phys. Rev. D 23 (1981) 1070. 
[12] A. Bodek et al, Phys. Rev. D 20 (1979) 1471. [26] 

[13] W.B. Atwood and G.B. West, Phys. Rev. D 7 (1973) 773. 
[14] T. Hamada and I.D. Johnston, Nucl. Phys. A 34 (1982) [27] 

382. 
[15] R.B. Wiringa, V.G.J. Stoks and R. Schiavilla, Phys. Rev. 

C 51 (1995) 38. [28] 

[16] O. Benhar and V.R. Pandharipande, Phys Rev. C 47 

(1993) 2218. [29] 



I. Niculescu, PhD Thesis, Hampton University, 1999. Un- 
published. 

I. Niculescu et al, Phys. Rev. Lett. 85 (2000) 1186. 
Y. Liang, M.E. Christy, R. Ent and C.E. Keppel (Jeffer- 
son Lab Hall C E94-110 Collaboration), nucl-ex/0410027 
http://www.jlab.org/~christy/cs_fits/cs_fits.html 
O. Lalakulich and E.A. Paschos, Phys. Rev. D 71 (2005) 
074003. 

O. Lalakulich, E.A. Paschos and G. Piranishvili, Nucl. 
Phys. B (Proc. Suppl.) 159 (2006) 133. 
Nuclear Physics B (Proc. Suppl.). 

O. Benhar, A. Fabrocini, S. Fantoni and I. Sick, 
Nucl. Phys. A579 (1994) 493. 

J. Carlson and R.Schiavilla, Rev. Mod. Phys. 70 
(1998) 743. 

J.E. Amaro, M.B. Barbaro, J. A. Caballero, T.W. 
Donnelly, A. Molinari and I. Sick, Phys. Rev. C 
71 (2005) 015501. 

T. Sato and T.-S.H. Lee, Phys. Rev. C 54 (1996) 
2660. 

T. Sato, B. Szczerbinska, K. Kubodera and T.- 
S.H. Lee, Nucl. Phys. B (Proc. Suppl.) 159 (2006) 
141. 

V. Tvaskis, J. Steinman and J. Bradford, Nucl. 
Phys. B (Proc. Suppl.) 159 (2006) 163. 

www.jlab.org/~kuhn/BoNuS_Welcome.html 



